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The energy route to the equation of state of hard-sphere fluids is ill-defined since the internal 
energy is just that of an ideal gas and thus it is independent of density. It is shown that this 
ambiguity can be avoided by considering a square-shoulder interaction and taking the limit of 
vanishing shoulder width. The resulting hard-sphere equation of state coincides exactly with the 
one obtained through the virial route. Therefore, the energy and virial routes to the equation of 
state of hard-sphere fluids can be considered as equivalent. 



PACS numbers: 



Given a fluid of particles interacting via a two-body po- 
tential (p(r), its equation of state (EOS) can be obtained 
in terms of the radial distribution function g(r;p,P), 
where p is the number density and (3 = 1/ksT is the 
inverse temperature, through a number of routesi The 
most common ones are the virial route, 

^ = Z(p,/3) = l+2«*- 1 « d p r drr d y{r-p^)^e-^\ 
P Jo or 
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the compressibility route 
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and the energy route 
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In Eqs. 0-©, P is the pressure, Z is the compress- 
ibility factor, d is the dimensionality of the system, 
Vd = (7r/4) d / 2 /r(l+c;/2) is the volume of a d-dimensional 
sphere of unit diameter, \ is the isothermal suscep- 
tibility, u is the internal energy per particle, y(r) = 
exp[(3ip(r)]g(r) is the cavity function, and h{r) = g{r) — 1 
is the total correlation function. If the exact function 
g(r; p, (3) is inserted, the three routes are thermodynam- 
ically consistent, i.e., 



X- 1 (p,(3) = ^- p [pZ(p, ( 
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On the other hand, if an approximate function g(r; p, (3) 
is used, the compressibility factor obtained directly from 
Eq. does not necessarily coincide with that obtained 
from the combination of Eqs. (J2J and (@J or with that 
obtained from the combination of Eqs. © and JSJ. 



Let us now particularize to the hard-sphere (HS) in- 
teraction 




In that case, y(r; p, (3) = yHs( r ; po~ d ) and Eq. Q becomes 



Zns(po- d ) = 1 + 2 d - 1 v d pa d ms {a-pa d ). 
On the other hand, Eq. © reduces to 
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Since Zhs is independent of temperature and i*hs is m " 
dependent of density, Eq. (JSJ is trivially satisfied as = 
and so it is not possible in principle to get the compress- 
ibility factor from the internal energy. 

A way of circumventing the ill-definition of the energy 
route to the EOS of an HS fluid consists of considering 
a convenient interaction potential that encompasses that 
of hard spheres in certain limits but for which the tem- 
perature plays a relevant role. The simplest choice for 
such a potential is perhaps the so-called square-shoulder 
(SS) potential^ 



Vss(r-) 
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where e is a positive constant. For this potential, Eq. @ 
becomes 
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w 
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The SS interaction is equivalent to an HS interaction of 
diameter er in the infinite-temperature limit (/3e — > 0) 
and to an HS interaction of diameter a 1 in the zero- 
temperature limit {(3e — > oo). Of course, the equivalence 
SS — > HS also holds in the limit of zero shoulder width 
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Let us now imagine that yss(i", P-, P) is known, either 
exactly or approximately (e.g., as obtained from the 
Percus-Yevick approximation), and thus it is possible 
to compute uss(p,P) from Eq. l|10fl . Then, we can get 
^ss (p, 0) from Eq. J5J as 
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where we have made use of the property 
lim^e^o Zss(p, P) = Ztis(pcr d )- Taking into account 
the complementary condition lim^^oo Zss(p, 0) — 
Z KS (p<j' d ), Eq. dm yields 



Zm(pv' d ) ~ Z m (pv d ) 2 d ~ 1 dv d e 8 
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Equation (|12|l can be considered as the condition defin- 
ing the compressibility factor of an HS fluid associated 
with the energy route. To proceed further, we take the 
limit a' — > a. Thus, 
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a' — >(7 



Z^{po' d ) - Z m (pa d ) 
pa' d - pa d 



<J- d -^-Z m (pa d ), (13) 
dp 



Finally, integrating over density and imposing the ideal 
gas boundary condition Zhs(O) = 1, the virial EOS 
is reobtained. 

The generalization to mixtures (either additive or non- 
additive) is straightforward. In that case, the SS poten- 
tial for each pair ij is defined by Eq. with the changes 
e — ► £y, (j — * (Tij , and a' — > e>y = A<7y , where the factor 
A is common to all the pairs. The equation equivalent to 
Eq. I|12|) in the case of mixtures is 
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where {xi} are mole fractions and a d s = - X i X j 7 7 

(j'itf = Y,i,j x i x j a 'ij = ^ d<7 cff- Taking the limit A — > 1 in 
Eq. l(ho)l one gets 

Zns{pv d fi) = 1 + 2 d - 1 v d pJ^x i x j af j y% s (a ij ;pa d s ), 

(17) 

which is not but the virial EOS for an HS mixture. 

In summary, in this note I have shown that the ill- 
definition of the energy route to the EOS of an HS fluid 
is saved by first considering an SS fluid and then tak- 
ing the limit of vanishing shoulder width. The resulting 
EOS coincides exactly with the one obtained through the 
virial route. From that point of view, the energy and 
virial routes to the EOS of HS fluids can be considered 
as equivalent. 



lim 



drr d 1 y ss (r;p,f3) = -y ns (a; pa d ). 



Therefore, Eq. (|12|l reduces in the limit a' — > a to 



(14) 



-^Z m (pa d ) = 2 d - 1 v d ^-pcr d y m (<j;pa d ). (15) 
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